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Recently in Ref. [l| has been proposed a GUT model for ferniion masses and mixings with 
spontaneously broken S3 discrete flavor symmetry, where S3 is the permutation group of 
Cm ■ three objects. The S3 breaking pattern in the quark sector is not studied and need further 

^SJ ' investigation. Since in such a model the number of free parameters is greater than the num- 

> ; 

\^^ I ber of experimental observables, an analytical fit of all the parameters is impossible. To go 



o 



X 



forward with the model building and to deal with this problem we have used a statistical 



'^ . analysis. We have found that S3 is totally broken and the up- type quarks matrix is approx- 

^O ' imatively diagonal while down-type quarks matrix is not symmetric and it is parametrized 

o ■ 

by three couplings, g'^, g'l and 53. It has been found that g^ is slightly smaller than g'^ = 1 
Qh! and it is of order one, while ^3 ~ A^ where A is the Cabibbo angle. An analytical study of 



Mh. the dependence of Vcb and Vub from the couplings 5^ and gf is also presented. 
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I. INTRODUCTION 

In third quantization J2l] we can explain the origin of the fermion families and the origin of 
the permutation symmetry between them. Third quantization is an extension of the quantum 
field theory where the fields are treated as well as particles in second quantization. Discrete 
flavor symmetries, like the group of permutation of three objects 6*3, are very interesting since 
they can naturally explain neutrino mixing angles and mass hierarchies, for instance see Ref. [3|]. 
Some recent reviews in the study of hierarchical structure of fermion masses in the Standard 

J] . Indications toward grand unification gauge theory 



Model and its extensions are reported in Ref 
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(GUT), is the phenomenological tendency to unify of the gauge couphngs, and the theoretical 



impUcit possibihty to explain charge quantization and anomaly cancellation. In Ref. [ll] it has been 
proposed a model based on Eq x S^ in order to arrange quark, lepton and neutrino masses and 
mixings. Minimal GUT models give strong relations between quarks and leptons Yukawa matrices. 
For example in minimal SO(IO) models up quarks Yukawa matrix is equal to the Dirac neutrino 
Yukawa matrix. However quarks and neutrinos mass and mixing hierarchies are very different. In 
fact, neutrino oscillation data show a maximal atmospheric angle, large solar angle and ^13 close 
to zero, see Ref. [^l- In contrast the greatest angle in the quark sector, namely the Cabibbo one, is 
smaller than the corresponding leptonic mixing angle, see Ref. [6|]. Phenomenological quark/lepton 
difference seems in contrast with GUT features. Even though in Eq models it is possible to make 
a distinction between charged fermions and neutrinos as showed in Ref. [l|. Authors evidenced 
this possibility by means of the two standard model singlets contained in Eq. In such a model 
quarks Yukawas only emerge after 6*3 symmetry is broken by not renormalizable dimension five 
operators, while renormalizable operators only contribute to Dirac neutrino Yukawas. However the 
study of 5*3 breaking pattern in the quark sector has not been accomplished in Ref. [l|] and it will 
be performed in this paper. 

We will show that in our model up and down-type quark mass matrices are parametrized both 
by five complex couplings and three real parameters. Since the number of free parameters in the 
quark sector of the model is larger than the number of experimental data, a direct fit is impossible. 
One possibility to deal with this problem, is comparing the theoretical expectations with the 
experimental data by means of a montecarlo method (see Ref. iP\) that assigns a probability 
weight to points in the parameter space of the model according to a goodness of fit criterion. 
In this way we find the statistically preferred textures of the mass matrices for up and down-type 
quarks. In particular we select a large number of points in the parameters space initially uniformed 
distributed. After applying the experimental constrains, not all the selected points agree with data. 
To apply the experimental constrains it is been used the x^ function defined as below 

X^(0^^)=E( ^e.r ) (1) 

where O^^^ and cr^^^ are respectively the experimental data and their standard deviations while 
Qth g^j,g |.]^g corresponding values predicted by the model. The points in the configuration space are 
statistically taken with a probability proportional to exp(— x^/2). After experimental constrains 
are applied we get a not uniform distribution in the parameters space. Regions very high populated 
will be interpreted statistically more probable than regions with low density population. This 



method selects models statistically preferred by data giving the magnitude of the fitting parameters. 

The paper is organized as follows: in section [TI] we give the up and down quark mass matrices 
textures derived by Eq x 5*3 model, in section HIH we introduce the numerical method used and we 
report the results, in section ITVl we give the conclusions. 

II. S3 SYMMETRY AND THE QUARK MASS MATRICES 

The discrete group of permutation of three objects, called 5*3, contains six elements and it has 
four irreducible representations, namely a doublet 2, one symmetric singlet 1 and one antisymmet- 
ric singlet 1', see for instance Ref. [8|]. The product of two doublets is 

2x2 = 2 + 1 + 1'. (2) 

In particular, if (xi5X2) ^^'^ (Xi)/t2) ^^'^ ^'^o ^3 doublets, their product is 

XD = -Mxl X? + X? X2 , X\ x\ - X2 xl) 

Xs = -7^{X1X\ + X2XI) 

XA = -^(X2 X\-X1X2) 

where XD is the doublet of eq.([2]) while xs and XA are respectively the symmetric and the anti- 
symmetric singlets. 

Let be ■(/'i; "02 and ^3 three fields belonging to a triplet reducible representation 3 of S3. Then 
the triplet reducible representation decompisition is 3=2+1: 

V'S = ^(V'l +^2 + -03) 

ijD = (4f (-^2 + V'S) , -^(-2 Vl + V'2 + V'3)) (3) 



where ips is a symmetric sing^ 

In Eq X S3 model of Ref. 
below 



et of S3 and TpD is a doublet of 53. 

l| the quarks Yukawa matrices come from dimension five operators 



Lvuk = Yl {^ijk HuQLiURj4>l + gfjk HdQLidRj^ij + h.c. (4) 

where i,j, k = 1,2, 3, the scalar fields (/>"' are Standard Model singlets, Hu and Hd are electroweak 
doublets, Qii are the left-handed fermion electroweak doublets while URi and dRi are right-handed 



fermion electroweak singlets and oj*', are 5*3 tensors. We assign the left-handed fermions, the 
right-handed fermions and the scalar (/>^ ' to triplet reducible representations of S^, namely 

Ql = {ui,di)L r^ 3 , URi r^ 3 , dRi r^ 3, 4il' ~ 3. 

The Higgs electroweak doublets are ^3 symmetric singlets, i/u,d ~ 1- With this assignment both 
the operators in eq. (HD transform with respect to the permutation group as follows 

3a X 3b X 3c X 1 = (2a + la) X (2b + lb) X (2e + Ic) X 1 = (5) 

= (2a X 2b X 2c + 2a X 2b X Ic + 2a X lb X 2c + +la X 2b X 2c + la X lb X Ic) X 1. 

First four terms of five in eq.dS]) contain the symmetric and the antisymmetric singlets, see eq. ^. 
Only symmetric singlets give S^ invariant interactions because we have chosen the Higgs fields in 
symmetric singlet representation and 1x1' = 1'. Consequently Yukawa interactions in eq. ^ give 
only five terms S^ invariant [I4] both in up and down quark sectors. 
The lagrangian is, see also Ref. [9|, 

Lyuk = (4 QDdD4>D + 4 QDds<t>D + 4 QsdD(t>D + 4 QDdD<t>s + 4 Qsds4>s) H'^ 

+ (C^ QDUD(t>D + C2 QdUS^D + 4 QsUD(t>D + 4 QdUd4>S + 4 QsUS(ps) ^" 

+ h.c. (6) 

where Qd, do, ud, (pD, Qs-, ds, 4>s and ug follow from relations in eq.([3|) and c^' are arbitrary 
complex constants. With the following redefinition of the complex constants 
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"^1 = Y 3 (5 - 52 - 5L - S-R + 53) 

C2 = yl (g- 92 + '^ gi- gR- gs) 

•^3 = yl (g- g2- gL + 2 gR- 53) (J) 

C4 = yl ig + 2 §2- gi- gR- gs) 

C5 = 75 5' + 2(52 + gL + gR + 53) 

and using relations in eq. ^, from eq. ^ we get the following down quarks lagrangian 

Y, [gdldRH^i + g2dld]iH4>j + gLdl4^H<Pi + gRdld^j^H<j)j + g^dl^j^H <j)4) + h.c. (8) 



and the corresponding mass matrix is 
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where v^ = ((/)f) and v^ = {H'^). We observe that M"^ is parametrized by five complex coupHng 
constants and three scalar vevs, namely 13 real parameters. Similar lagrangian and mass matrix 
can be obtained for up quarks. 



We want to fit the hierarchies of the coupling constants §2, gi, 9r and 53 13], then we 
parametrize their absolute values as power of the Cabibbo angle as below 



u,d u,d \t"' 

92 = ^2 ^' 



9l 



at'' A*r 



u,d 

9r 



ai'" x^r 
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ag"'" ^3 



(10) 



where A = 0.2, the t exponents are free real parameters and a^' , a^' , 04' and a^' are complex 
phases. 



III. THE NUMERICAL FIT 



The method 



The goal is to extract the values of tg' > ^r ; ^^ and t^' of eq. (jlOp . namely the magnitude 



order of the gs couplings, and the values of the vevs of the scalars v^ 



u,d 



and V3' from the 



experimental measurements. A direct fit of the data is not possible, since the number of free 
parameters in eq. Q is much larger than the number of observables, six mass eigenvalues plus 
four CKM parameters. The main obstacle comes from the coefficients a"' , whose phases are not 
theoretically known. To cope with them, we will treat this uncertainty as a theoretical systematic 
error. Namely, we have assigned a flat probability to all the coefficients a"' with 



|a-' I = 1 , < arg(aj' ) < 27r. 



(11) 



The exponents and the vevs are chosen in the ranges 



,u^d ,u,d 
^2 ' ''/J ' 



,u,d 



,u,d 

h 



U<Io , tjy , tj , to < 



ry , u,d u,d , u,d -, 

<v^' , V2' < V3' = 1, 



(12) 



and we randomly take them with a flat distribution in logarithmic scale, v^' , V2' and v^' must 
satisfy the above constraints since (by definition) we choose the entry (3,3) of the matrices in ([9]) 
to be the largest one. For any random choice of the coefficients a^' , of the exponents t and of 
the vev's v'^' we respectively get two numerical matrices for the up and the down sectors. The 
diagonalization of these two matrices, gives six eigenvalues, corresponding to the physical up and 
down quark masses to be compared with the experimental values of the masses runned at the 
unification scale 2 • W^^Gev (see Ref. [lO"] ) and reported in Table HI The multiplication of the two 

Mass Reference value 

m,(MeV) 0.835lt°:i636 

TOc(MeV) 2A2M76tllfolt 
mt(GeV) 75AM8tlil^, 
md(MeV) 1.7372t;;:^^^^ 
TO^MeV) 34.597lt5:i9n 
ma(GeV) 0-957Atroml 

TABLE I: Quark masses runned at 2 • lO^^Gev scale in nonSUSY standard model. 

unitary matrices that diagonalize on the left the up and down mass matrices respectively, yields 
the quark mixing CKM matrix. The CKM is parametrized by the four Wolfenstein parameters A, 
A, p and r] and their experimental values (see Ref. [g]) are reported in Table HIl We have collected 



CKM parameter 


Reference value 


A 


0.2237 ±0.0033 


\Vcb\ 


(41.0 ±1.6) X 10-3 


P 


0.225 ±0.038 


V 


0.317 ±0.041 



TABLE IL CKM parameters. 

a large statistical sample of events. Each one event can be compared with the experimental data 
through the x^, namely the event is accepted with probability 

P(Of)=e-5X^(or) (13) 

where the x^ is defined in eq. ([T|). Before applying the experimental constraints, events are ho- 
mogeneously distributed in the variables t"' , f"' , and probability distributions are fiat, but after 
applying the weight corresponding to eq. (Iisp . only points lying in well defined regions of the space 



t"' and v^' have a good chance to survive and the events are not uniformally distributed in the 
configuaration space. The density is related to the probabihty that an event predicts masses and 
mixings in a range compatible with the experiments. Even if our montecarlo approach gives more 
predictive and accurate models, we also emphasize that one should not think these results as true 
experimental measurements. They only give "natural" range of values for the exponents f"' and 

uA 



B. The result 



The statistical analysis gives different regions of points in the parameters space accordingly with 
the distribution exp(— x^/2). Due to a numerical instability we are unable to compare different 
regions of solutions. Whatever isolated points in the configuration space are statistically disfavored 
since unphysical phases in eq. (jlOp should be tuned in order to fit the experimental data. While 



points in regions with high density population are statistically preferred by data as explained in 
the previous section and in Ref. 7|. 

Let us now to consider the scalar vevs obtained with the montecarlo analysis. In Fig. [1] are 
plotted on the left up-type quark vevs and on the right the down-type quark vevs, in particular it 
is plotted first generation in the y axis vs second generation in the x axis. The third generation of 
vevs are fixed to be Vo = 1 and vi = 1. From Fig. [1] we have the bounds below 



1 1e-05 




E 0.001 



1e-04 0.001 

second family vevs 




0.01 0.1 

second family vevs 



FIG. 1: Values of the VEV of the scalar SU{2) singlet fields in the up on the left and down sectors on the 
right. 



10"" < v\ < 10"* 
2 • 10"^ < v^ < 3 • 10"^ 



10"^ < vf < 10^2 
2 • 10^2 < yd < 6 . 10- 



in agreement with the observed mass hierarchies at the unification scale (see Table U]) 



Vi : V2 : v^ = niu ■■ rUc ■■ mt = X^ : X^ : 1, 
vf : V2 ■■ v^ = nid : rus ■■ rub = X'^ : X^ ■■ 1. 



Relations ^^ and 031) show that in quark sector the S3 symmetry is totally broken. 



(14) 
(15) 



The magnitude of §2' , g^' , g^ and gfg' are related to the ts exponents like in eq. (jlOp . 
In Fig. ([2]) we report the montecarlo results for the ts exponents given by the statistical analysis 
allowed after the application of the experimental constrains eq. (J13p . 
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FIG. 2: i exponents that parametrize the magnitude of the g couplings in the general case. We have plotted 
the plane i4,t%), {tl,t\), {^,4) and {t^,4). 



Prom these figures it is not difficult to obtain the following bounds for the t exponents 



0< t^ < 
4.5 < q < 



3 < t^ < 8, 
5 < q < 8. 



(16) 



< t^ < 1 ; 2.5 < 4 < 8, 



(17) 

3^o, u^ (-2 "° 



2.5 < f^ < 3 ; 5 < 4 < 



The upper bounds t =8 in eq. (|T6|) and eq. (|T7j) is the cutoff used in the montecarlo method, see 
eq. (|12p . An upper bound t < a < 8 correspond to a g coupling lower bound, namely \g\ > A". In 
the up quark sector, gs couplings are not bounded from below since the corresponding ts exponents 
could be large, namely in eq. (J16p all the ts exponents can be taken equal to the cutoff used in the 



statistical analysis. This means that all the gs couplings in the up quark sector can be very small 
of order A®. In contrast, in the down sector the couplings g^ and g^ are bounded from below[lJ] 

\gi\>X, |53l>A^ (18) 

while the couplings (72 ^^^'^ Sr i^ the down sector can be very small like in the up quark sector. 
Experimental data statistically prefer models where the absolute value of the |(7^| coupling is large, 
see eq. ()18p . In particular |(/^| can be close to one or to the maximal value allowed in the montecarlo 
analysis, see eq. (|12p and it is slightly smaller than g'^ = 1. Indeed from relations (J17p and (J18p 
the gf coupling also has an upper bound, that is \g'^\ < A^'^ ^ A. Therefore from the statistical 
analysis we get the following hierarchies between the gs couplings 

J- — 5 —9 — \9l\ ^ \9'i\ ^ \92h \9Rh \92\1 \9Lh \9Rh 193 I- \^^) 

In Table IIIII we report a numerical example of values for the gs coupling and for the vevs that fit 
good data. 



up down 



0«/0| = TO„/mt </(/)f = 0.000077 
(\)ll<i)l = mc/mt (pi/ipi = 0.054783 
3 = 1 g = 0.848 

51,-0 gL = 0.305 e"* ^-^"^ 

gn ~ gR = 0.002 e"* o^^ 

gL^O gL ^ 0.009 e' ^-^^ 



TABLE III: Example of set of values of g couplings and vevs that fit the data. 
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C. The CKM and the couplings gf and g^ 



The statistical analysis results in previous section, show that in first approximation the experi- 
mental data can be fitted only with couplings g^, g , gj^ and gf couplings (5"' are fixed), besides 
the vevs of the scalar fields. In this section we want to study the relation between the gf^ and g'f 
and the entries of the CKM mixing matrix. Accordingly we fix the values of the scalar vevs like 
in eq. (fT4|) and eq. (fTSJ) . Prom eq. (fT4|) . eq. (fTSJ) and eq. (fT9]) we get, up to correction of order A^, 
the mass matrices below 

^ A^ giX^' + gi giX^' + giy?^ 
^H gi\' + 9i A2 gi\^ + gi\' ■ (20) 
\gi+gi\^ gi+gi\' 1 / 

The up quark mass matrix in eq. (j20l) is almost diagonal and the quark mixing matrix Vckm 
is approximatively given by the unitary matrix Ud that diagonalize on the left the down quark 
mass matrix eq. (I20p . namely f/j M^ Mj Ud = D^. We observe that in general the g'l and gf are 
arbitrary complex variables. In order to simplify the problem, we assume g'l and gf reals and the 
CKM complex phase r] is not fitted. 

In the standard parametrization Vckm = R23 ^13 R12 where Rij are rotations in the {ij) plane 



M" 






( X' o\ 

A^ 
Vo ij 



of 9ij angles, then \Vcb\ ~ sin ^23 and \Vub\ ^ sin^is. From mass matrices (I20p we have 



IT/, I ~ 



\ 



1 + 



il + 2g}W-l) 



(21) 



'l6ff|A4 + (A4 - 1)2 + 4^4 (1 _ A4 + AS) 
where gL = A*^, so in first approximation \Vch\ is only a function of the coupling gf^. In Fig. [3] we 
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FIG. 3: Plot of the function \Vcb\ = K&I^l) and |K&| = Vub{h) 



plot the function \Vcb\ = Vcb{tL)- We observe that when ti — 0.5 the corresponding \Vcb\ agree 
with data and ti — 0.5 is close to the value from the montecarlo, see in eq. pT|l . 
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Analogously we can use the remaining free parameter g^ to fit the \Vub\ as below 



V2V V9 + 53" + 2453 A2 - 8<7i + 253^(10A4 + 4A2 - 1) " ^ ^ 

In Fig. [3] we have plotted the function Vubits)- l^fel agrees with experimental data when ^3 ~ 3, 
in agreement with statistical analysis. \Vub\ is approximatively only a function of g'l coupling. 
Assuming for instance t^ = 0.55 and t3 = 3.1 in eq. (pOj) . we have \Vus\ — 0.21, \Vub\ — 0.004 and 
\Vcb\ — 0.043 that agree with data. 

IV. CONCLUSIONS 

Recently it has been proposed in Ref. [l| a model based on Eq x 5*3 in order to explain fermions 
mass and mixing hierarchies. However a full study of 5*3 breaking pattern was not accomplished. 
We carry out such a study. In the model of Ref. [l|] quark mass terms appear as higher order 
dimension five operators after 5*3 is broken. We have showed that in such a model quark mass 
matrices are parametrized by five complex couplings 5"''^, §2' , g^' , g^ , 53' and three real 
parameters v^' , V2' , v^' that are the vevs of the scalar fields that break S3. Even if we can 
assume g^ = g = 1 and t"" = v^ = 1, the number of free parameters in the model is larger than 
the number of experimental data then it is not possible to fit the parameters. To go forward we 
have used a numerical approach finding the statistically preferred textures of the mass matrices 
according to a goodness of fit criterion. We have found that the permutation symmetry 5*3 is totally 
broken and vev hierarchies should be v^' <C ^2' *C ^3' , and data statistically prefer solutions 
with large {g'll coupling, namely slightly smaller than g"^ = 1, and 153! ~ A'^ while the remain 
couplings in up and down quark textures, are negligible and are of order A^. The resulting up 
quark mass matrix is almost diagonal and the quarks mixing matrix is approximatively given by 
the down quark mass matrix that is not symmetric. The hierarchies between the yukawa couplings 
is 

1 = 5" = /=^ Igil » \9i\ » bfl, \9rI \92l \9l\, KrI m- (23) 

Ultimately we have studied the relation between the couplings g'j^ , g^ and the CKM mixing matrix. 
We have found that Vcb is approximatively a function of gf coupling and Vub is approximatively a 
function of g^ coupling. 

The study of the breaking pattern of the S3 symmetry in the quark sector, gives indications to 



go forward in the model building. In fact recently in Ref. 



llj l it has been studied a model that could 
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explains quark, lepton and neutrino mixings and masses, making use of the montecarlo statistical 
analysis results reported in this paper. In particular in the model of Ref. 



Ill ] are explained the 



hierarchies in eq. (123p between the gs Yukawa couplings. In third quantization models it is possible 
to extend the concept of family to the gauge bosons through the semidirect product of a gauge 



group G, with the group of the permutation of n objects, namely G" >iSn- In Ref. [ll| it has been 
proposed a Eq X S4 model where Eq is the grand unified gauge group and 54 is the permutation 
symmetry of four objects. It is showed that embedding 53 C 5*4 it is possible to better explain 
neutrino oscillation, while Eq could be useful in order to understand the hierarchies between the gs 
couplings reported in eq. (|23]) . It is presented a possible breaking pattern of Eq >i S4 in agreement 
with our montecarlo results. In the following we report such a breaking pattern. 

Assume that the starting group Eq yjS^ is broken into the group Gi = {SM x U{l)r x C/(l)t)^ > 

< S'4 where U{l)r and U{l)t are defined as Eq D 50(10) x U{l)t D SU{5) x U{l)r x U{l)t and 

SM is the Standard Model gauge group. At some scale Mi the group Gi breaks into the group 

G2 = {SM X U{l)r X U{l)t)'^ X S3. It is easy to see that the only operators compatible with G2 

are 

The coupling g will be of order of magnitude g ~ Mi /A where A is the characteristic scale of 
Sg XI 54. At another scale M2 < Mi the group G2 breaks into G3 = {SM x ^(l)down)^ >< 5*3 
where C/(l)down is defined as the linear combination of U{l)r and U{l)t so that right-handed down 
quarks do not carry f/(l)down charge. At this scale also the operators 

Y,9Ldld'nH'ct>i 

hi 

are allowed and gL — M2/A. If Mi ~ M2 it is possible to explain the relation g c^ gi in eq. (j23p . 
At the scale M2 the operators proportional to ^2; di ^^^ Br ^^^ ^ot allowed. They only appear 
after the breaking of C/(l)down at some scale M3 < M2. If M2 ^ M3 we can explain the hierarchy 

gf » gi gi, gi 
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